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Production and Evolution of Perturbations of Sterile Neutrino Dark Matter 



O 
O 

(N 



> 
O 
m 



in 
o 



6 ■ 



X 



Kevork Abazajian 

Theoretical Division, MS B285, Los Alamos National Laboratory, Los Alamos, NM 87545 

Sterile neutrinos, fermions with no standard model couplings [SU(2) singlets], are predicted by 
most extensions of the standard model, and may be the dark matter. I describe the nonthermal 
production and linear perturbation evolution in the early universe of this dark matter candidate. I 
calculate production of sterile neutrino dark matter including effects of Friedmann dynamics dictated 
by the quark-hadron transition and particle population, the alteration of finite temperature effective 
mass of active neutrinos due to the presence of thermal leptons, and heating of the coupled species 
due to the disappearance of degrees of freedom in the plasma. These effects leave the sterile neutrinos 
with a non-trivial momentum distribution. I also calculate the evolution of sterile neutrino density 
perturbations in the early universe through the linear regime and provide a fitting function form 
for the transfer function describing the suppression of small scale fluctuations for this warm dark 
matter candidate. The results presented here differ quantitatively from previous work due to the 
inclusion here of the relevant physical effects during the production epoch. 

PACS numbers: 95.35.+d,14.60.Pq,14.60.St,98.65.-r 



I. INTRODUCTION 

The nature of dark matter remains one of the most 
significant unsolved problems in cosmology and particle 
physics. The abundance of dark matter has been pre- 
cisely determined by observations of anisotropics in the 
cosmic microwave background 1] , and the measurements 
of the growth of cosmological structure in the clustering 
of galaxies 0] and in the Lyman-a forest 0. The fun- 
damental nature of the dark matter, however, remains 
unknown. 

One natural candidate is a fermion that has no stan- 
dard model interactions other than a coupling to the 
standard neutrinos through their mass generation mech- 
anism 0, . Due to their lack of interactions and associ- 
ation with the neutrino sector, such fermions are referred 
to as sterile neutrinos. Observations are consistent with 
sterile neutrinos as the dark matter for a narrow mass 
range for the standard production mechanism. In this 
allowed range of masses, the sterile neutrino has a non- 
negligible thermal velocity component, and is therefore a 
warm dark matter (WDM) candidate. 

The prevalent ansatz of an absolute cold dark mat- 
ter (CDM) component in galaxy formation is not strictly 
valid even for one of the most cited CDM candidate, the 
lightest supersymmetric particle, which has a small but 
non-zero velocity dispersion 0, Q- The damping scale 
at which thermal velocities of the dark matter cut off 
the growth of gravitationally bound structures remains 
an open question. One principal challenge to the CDM 
paradigm is the order of magnitude over-prediction of 
the observed satellites in galaxy-sized halos such as the 
Milky Way H, 0, E3, EH 

Warm dark matter suppresses 
dwarf galaxy formation, which may occur through frag- 
mentation of larger structures [T^ . Semi-analytic galaxy 
formation modeling has found that the number of dwarf 



galaxies formed in satellite halos may be suppressed due 
to the reionization, stellar feedback within halos, and/or 
tidal stripping of satellites [H Q EE E3 Such 
semi-analytic modeling is both powerful and malleable, 
and must be verified in robust hydrodynamic simulations 
of galaxy formation. Whether a minor or major suppres- 
sion of small mass halos is beneficial or detrimental to the 
suppression of dwarf galaxy formation remains unsolved. 

Four more problems in the CDM paradigm may ben- 
efit from the reduction of power on small scales from 
WDM. First is the reduction of the prevalence of halos in 
low-density voids in N-body simulations of CDM struc- 
ture formation, consistent the apparent dearth of mas- 
sive galaxies within voids in local galaxy surveys |12lll 8| . 
The second is the relatively low concentrations of galax- 
ies observed in rotation curves compared to what is pre- 
dicted from the ACDM power spectrum 0, [2(J , which 
can be relieved by a reduction of the initial power spec- 
trum of density fluctuations at small scales |2lL . The 
third is the "angular- momentum" problem of CDM halos, 
where gas cools at very early times into small mass halos 
and leads to massive low-angular momentum gas cores in 
galaxies, which can be alleviated by the hindrance of gas 
collapse and angular momentum loss through the delay of 
small halo formation in a WDM scenario [23| . The fourth 
problem is the formation of disk-dominated or pure-disk 
galaxies in CDM models, which is impeded by bulge for- 
mation due to the high merger accretion rate history in 
CDM models, but may be alleviated with WDM pll^ . 

Sterile neutrinos arise naturally in most extensions to 
the standard model of particle physics. Singlet neu- 
trinos with masses relevant to oscillation experiments 
and dark matter can arise in grand unified models |'26| . 
string- inspired models |2?|. and models with large ex- 
tra dimensions psl l29j . In reverence to Occam's ra- 
zor, the minimalist model of Ref. |30j produces the 
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neutrino-oscillation inferred neutrino mass hierarchy, the 
dark matter via a sterile neutrino, as well the observed 
baryon asymmetry in a neutrino minimal standard model 
(i/MSM) that introduces Majorana and Dirac neutrino 
mass terms to the standard model Lagrangian. There 
is also an indication that one or more light sterile mass 
eigenstates may cause the flavor transformation seen in 
the Los Alamos Liquid Scintillator Neutrino Detector 
(LSND) experiment 0, Another intriguing moti- 
vation for the presence of a dark matter sterile neutrino 
is abundance of anomalously high pulsar velocities that 
may be difficult to produce in the convective hydrodyam- 
ics in a supernova, but may be produced in assymet- 
ric sterile neutrino emission from a hot nascent neutron 
star [s^L IsEI]. However, whether convective overturn 
and the resulting global asymmetry in the ejecta alone 
can power the observed distribution of pulsar velocities 
remains an open question |36l l37j . 

The potentially beneficial effects of the suppression 
of cosmological small-scale structure in WDM can also 
lead to observational conflicts if the suppression extends 
to excessively high mass and length scales. As I shall 
show in detail below, the suppression scale monotoni- 
cally decreases with increasing sterile neutrino particle 
mass. The reionization of the universe by a redshift of 
z ~ 6 requires sufficient structure formation at very early 
times, and can place one of the most st ring ent lower- 
bounds on the sterile neutrino mass [38l l39( . The ra- 
diative decay of the sterile neutrino dark matter may, 
however, increase the hydrogen ionization fraction, aug- 
menting molecular hydrogen formation , ga s cooling, star 
formation, and therefore reionization [4C|. One of the 
best direct measures of clustering at small scales is the 
clustering observed in intervening gas along the line-of- 
sight to a quasar, known as the Lyman-a forest [4l| . 
Statistically-consistent constraints allowing freedom in 
all cosmological parameters and constraints from the cos- 
mic microwave background, galaxy clustering, and a mea- 
surement of clustering in the Lyman-a forest gives a lower 
limit for the sterile neutrino dark matter particle mass as 
m s > 1.7 keV (95% CL) H2- 

In all studies of the types discussed above of galaxy 
formation, cosmological reionization and clustering in the 
Lyman-a forest, robust conclusions require an accurate 
initial description of the fluctuations arising from the 
early universe. In Jill I outline the production mecha- 
nism of sterile neutrino dark matter through the varying 
particle population in the early universe and the QCD 
transition, and describe the resulting non-thermal sterile 



neutrino energy distribution. In TOII I follow the evo- 
lution of the non-thermal sterile neutrino dark matter 
perturbations through the radiation dominated era into 
the linear regime of the matter dominated era. This 
is related to the CDM perturbation spectrum through 
a transfer function. My results differ significantly from 
previous work on sterile neutrino dark matter perturba- 
tion evolution |43J, |4J] which neglected the effects of the 
changing particle population in the early universe, the 
QCD transition, and the dilution of the dark matter due 
to annihilation. 

The sterile neutrino particle dark matter candidate 
studied here may be embedded in several extensions be- 
yond the standard model of particle physics, as a super- 
partner, or it may have properties that would have it 
couple to other species at higher temperatures, including 
the inflaton. However, knowing its behavior in the early 
universe in such extensions would involve introducing a 
much more model dependent interaction and production 
mechanism, as well as knowledge of all of the degrees of 
freedom that are present and may annihilate and dilute 
the sterile neutrino in such a model between the coupling 
epoch and today. Even in such models with higher en- 
ergy scale couplings, dilution may render the abundance 
of sterile neutrinos negligible upon entering the produc- 
tion epoch considered here of T < 300 MeV. The pro- 
duction mechanism studied here is a minimalist extension 
to the standard model through the neutrino mass gener- 
ation mechanism, and does not require what would be at 
this point speculation of higher energy physics. There- 
fore, the initial abundance of sterile neutrinos entering 
the oscillation-production epoch is taken here to be nil. 



II. PRODUCTION 

Sterile neutrinos of interest for dark matter are never 
coupled to the primordial plasma. The production of 
sterile neutrinos in the early universe within the mass 
range of interest for warm to cold dark matter occurs at 
temperatures where collisions dominate the evolution of 
the neutrino system, and matter-effected oscillations are 
suppressed by induced thermal masses and the quantum 
zeno effect. The true time evolution of the system is de- 
scribed by that of the density matrix [45|. However, the 
collision dominated regime allows a simplification of the 
density matrix evolution to a quasi-classical Boltzmann 
equation of the form 



d f( „ d 1 r a (p)A 2 (p)sin 2 20 
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with a corresponding evolution equation for the antineu- tive and sterile neutrino distribution functions, as a func- 
trino distributions. Here, f s (p,t) and f a (p,t) are the ac- 
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tion of momentum, p, and time, t; H = a/a is the Hubble 
expansion for scale factor a; A(p) = Sm 2 /2p is the vac- 
uum oscillation factor dependent on the mass-squared 



difference Sm — 



between the active and sterile 



neutrinos; the mixing angle between the two flavor states 

is e. 

The production is driven by the collision rate 



r tt (p) 



1.27 GlpT 4 , 
0.92 G|pT 4 , 



A 4 ) T ! 



(2) 



and is augmented at temperatures, T, above the quark- 
hadron (QCD) transition, where quarks and the massive 
leptons fi and r contribute [46j, and are included in the 
calculation here. The system is damped by the quantum 
Zeno effect, D{p) — T n (p)/2, and mixing is suppressed 
by the thermal potential 47] 



V T ( P ) = 
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which has contributions from thermally populated lep- 
tons of the same flavor as the active neutrino. 
The asymmetric lepton potential for flavor a is 



V L = V2G f 



2 (n Va + (n Vp - n Vf3 ) - 



(4) 

Here, I will only consider lepton number symmetric uni- 
verses; however, the limits on the lepton asymmetry al- 
low for a non-negligible lepton number |48t |49j, |50| which 
can drive resonant sterile neutrino dark matter produc- 
tion pi, [1^ . One can neglect the asymmetric potential 
due to the baryon number n n , which is subdominant at 
all temperature scales of interest for sterile neutrino dark 
matter production. The temperature of peak production 
is approximately 



(1 ft \ 
3kev) 



1/3 



(5) 



The time evolution of the temperature-dependent ther- 
mal potential and collision rate needed to integrate 
the production Boltzmann equation requires knowing 
the time-temperature relation, and therefore the evolu- 
tion of the background plasma. Details of the general 
time-temperature relation are given in the appendix of 
Ref. [4r| . and are summarized here. The time tempera- 
ture relation depends on the expansion-dependent change 
of the temperature 



da _ dptot , v-i 

-jj, — ^rp IPtot +Pt0tJ 
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and the Friedmann equation governing the expansion 
rate, da/ dt = 3H . The standard evolution of the pressure 



and temperature are affected by the changing thermal 
population of particle species, which 1 calculate using the 
known standard model particle mass distribution 51]. 

A critical consideration for dark matter production is 
the fact that production occurs near the temperature 
of the quark- hadron transition 52|. Lattice QCD cal- 
culations of two massless u, d quarks and an infinitely 
massive s-quark show a cross-over type transition at 
Tqcd = 173 ± 8 MeV 53]; more realistic (2+1) quark 
flavor lattice calculations find a transition at Tqcd = 
169±12 (stat.)±4(sys.) MeV 54]. Overall, due to hadron 
and lepton population as well as the quark-hadron tran- 
sition, the statistical degrees of freedom of the plasma, 
(7* , change by nearly an order of magnitude in the tem- 
perature range of interest for sterile neutrino dark matter 
production. In addition, the annihilation of these species 
heats the coupled plasma relative to the (constantly) de- 
coupled sterile neutrinos, which are subsequently diluted 
and spectrally distorted. Note that production is also 
affected by the temperature range of the softness of the 
cross-over transition. Motivated by Ref. I model 

the softness of the crossover as rapid, over 5 MeV at 
Tqcd — 170 MeV. Changes to the softness of the transi- 
tion affect the production abundance by only a few per- 
cent. 

For low-mass sterile neutrinos, m s < 0.5 keV, pro- 
duction occurs at low temperatures where the statistical 
degrees of freedom in the primordial plasma are chang- 
ing minimally, and one may set to be constant so that 
the production Eq. JIJ can be integrated analytically 3 . 
In this approximation, the momentum dependence of the 
production disappears, and the sterile neutrino momen- 
tum distribution is simply a suppressed active neutrino 
distribution. 

However, for m s > 0.5 keV, in order to accurately in- 
clude the strong changes in the plasma background, the 
quark-hadron transition, dark matter dilution, the modi- 
fication of scattering rates, and flavor-dependent modifi- 
cation of the thermal potential [Eq. J3J)] , the production 
equation must be integrated numerically, and the result 
is momentum-dependent. In order to test the accuracy of 
the numerical calculation, I solve the analytically solvable 
case of a simple power-law time-temperature relation as 
well as all of the other simplifications of the analytic case, 
but with the full momentum distribution of neutrinos. 
The numerical integration recovers the analytic result 
within 1%. In Fig.^ I show the resulting relative distri- 
bution of sterile to active neutrinos p(e) = / s (e)// a (e) 
for masses in the range 0.3 < m s < 140 keV, with 
e = p/T. Note that all cases are distorted from a pure 
constant suppression, and are increasingly distorted for 
more massive neutrinos, which are produced at higher 
temperatures [Eq. JSJ] where the effects described above 
are more pronounced. One of the dominant effects is the 
cooling and subsequent enhancement of the low momen- 
tum sterile neutrino distribution due to the heating of 
the plasma, including the active neutrinos, from parti- 
cle annihilations. Another feature is the enhancement at 
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FIG. 1: The resulting relative distribution coming from the 
production epoch for sterile neutrino dark matter relative to 
the active neutrinos, p(e) — / s (e)// a (e), (e = p/T) over a 
mass range 0.3 < m s < 140 keV, for fifty cases, and increasing 
rris having decreasing distribution amplitudes. All cases have 
SIdm = 0.26. The upper thick (red) line is for the case of 
m s = 1.7 keV and lower thick (blue) line is for the case m s — 
8.2 keV. 



e < 1 for lower mass neutrinos due to the production of 
these momenta during the slow temporal evolution of the 
temperature through the quark-hadron transition. 

The work on sterile neutrino dark matter production 
of Ref. [55| follows that of Ref. 4j except for an exten- 
sion to m s > 0.5 keV via the inclusion of an unknown 
factor (g£/10.75) within the production "prediction" re- 
lationship, with an effective production statistical degree 
of freedom, g'^. The resulting relation lacked predictiv- 
ity due to the power-law dependence on an unknown . 
The major effect of the high degrees of freedom at the 
production temperatures is dilution of the sterile neutri- 
nos, so that the production dependence on an increased 

is an inverse relation. The sign and value of the de- 
pendency on gl in the production prediction equation of 
Ref. [55j was noted in Ref. |44( as a typographical error. 
Using the common but incorrect choice of gl — 10.75, the 
production relation in Ref. |55| is inaccurate, as shown 
in Fig. 

In Fig. [21 I show the contours of critical density in 
sterile neutrino dark matter density for varying mass 
and mixing for an electron neutrino flavor mixed with 
the sterile neutrino. I also show the results of Ref. [5{| 
for the common choice of g'% = 10.75, which is inaccu- 
rate, but choices of a more realistic effective production 
temperature g' t increases the discrepancy. I also show 
the results of Ref. 46], which used Tqcd = 100 MeV 
and a different model of the hadron and lepton popula- 
tion distribution based on an older catalogue ofparticle 
masses. The critical density contour of Ref. [46( lies at 
higher mixing angles because an increased coupling was 
required to offset the dilution of the QCD transition set 



FIG. 2: Contours of predicted critical density SIdm = 0.26 
from the direct numerical calculation in this work (solid 
red), the fit provided here (dot-dashed purple), the results 
of Ref. fltl which used a different quark-hadron transition 
and particle population model (dashed blue), and that of 
Ref. |55j| (dotted black) for the inaccurate but common choice 
of g' t — 10.75 (more realistic choices of g' t make the predicted 
abundance more inaccurate). Also shown are the upper flux 
constraint from X-ray observations of the Virgo cluster [EiJ, 
the constraint from the diffuse X-ray background [5^ . the 
lower mass constraint from the CMB, SDSS galaxy cluster- 
ing and Lyman-a forest, and the possible constraint from also 
including the high-resolution (HR) Lyman-a forest |4^ . 



at Tqcd = 100 MeV , which occurs below the bulk of 
production of all sterile neutrino masses m s > 1 keV. 
The full numerical results for the predicted dark matter 
abundance from this work are fit well by the relation 
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3.40 keV 



sin 2 29 



10- 



-0.615 



0.26 



x <^ 0.527 erfc 



1.15 



( T < 



QCD 



V170 MeV 



0.5 



2.15 



,(7) 



for an electron neutrino flavor mixed with the sterile 
neutrino. This expression is valid for 135 < Tqcd < 
300 MeV. Note that the quantity within curled brackets 
is unity for Tqcd = 170 MeV. 

This new relationship between the critical density frac- 
tion, m s and sin 2 29 modifies the flux constraint observed 
from the Virgo cluster 56], so that the inferred upper 
bound on the mass of the sterile neutrino dark matter is 
now 



< 8.2 keV. 



(8) 



This is the result of the fact that the constraint in 
Ref. [5(| is not a direct mass constraint but a flux 
constraint, which is related to the radiative decay rate 
through the mass-mixing angle relation. Lines of con- 
stant flux follow the fourth power of the mass since the 
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decay rate increases as the fifth power, but the number 
density in the field of view decreases proportionally with 
the mass. Using the production relationship, Eq. (JZJ), the 
diffuse X-ray background limit of Ref. [57j is 

/0 x 0.538 

m s < 8.89 keV (^) , (9) 

for central values of the cosmological parameters, and is 
shown in Fig. [21 The constraints from unresolved X-ray 
sources derived by Mapelli & Ferrara |5S| are similar to 
Eq. ©, when using the production relation Eq. Q) 59] . 

III. PERTURBATION EVOLUTION 

The standard cosmological model of structure forma- 
tion from adiabatic Gaussian fluctuations seeded by an 
inflationary epoch is affected by perturbation growth in 
the radiation through matter dominated eras. The distri- 
bution of velocities of the dark matter suppresses fluctua- 
tions below its free streaming scale, which increases with 
the mean dark matter velocities and decreases with its 
mass. Since sterile neutrinos are produced non-thermally, 
their full energy distribution must be included in an ac- 
curate calculation of the fluctuation spectrum arising 
from the linear growth epoch. I use the approach of 
the covariant multipolc perturbation evolution equations 
for massive neutrinos in Ref. |60j and implemented in 
the Code for Anisotropies in the Microwave Background 
(CAMB) The multipole equations depend on the 

value of the massive neutrino energy distribution and 
its momentum derivative, but I will not reproduce them 
here. 

I calculate the growth of perturbations through the 
radiation and matter dominated epochs of sterile neu- 
trino dark matter with CAMB. I include directly 
the numerically-calculated momentum-dependent sterile 
neutrino distribution functions and their derivatives from 
the solution of the quasi-classical Boltzmann Eq. as 
described in the previous section. The resulting linear 
matter power spectra today at redshift zero are shown in 
Fig. [21 for a range of sterile neutrino masses from 0.3 to 
140 keV, along with the related CDM case. 

A useful form of the suppressed perturbation power 
spectrum P ste riie(&) relative to the CDM case is a sterile 
neutrino transfer function of the form 

which can be used to convert any CDM transfer function 
to that of sterile neutrino dark matter. 1 find a fitting 
function that describes the transfer function of the form 

T s (k) = [l + (ak)T", (11) 

where 

— (l^v) S (^)W''- lM - (12) 
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FIG. 3: Shown are the resulting linear matter power spectra 
for nonthermal sterile neutrinos in the mass range 0.3 < m s < 
140 keV (gray/cyan). The thick (red) low-fc suppression case 
is for the lower-mass limit inferred from the Lyman-a forest 
(rn s > 1.7 keV) , and the thick (blue) high-fc suppression case 
is for the upper-mass limit from X-ray observations of the 
Virgo cluster (m s < 8.2 keV). The CDM case is the dashed 
(black) line. Measures of large-scale structure in the linear 
regime are in the region of O.Ol/i Mpc -1 < k < 0.2h Mpc -1 
for galaxy surveys, while neutral gas clustering observed in 
the Ly-a forest may extend observations of linear structure 
to O.lh Mpc -1 < k < 3h Mpc" 1 . 



and a = 0.189, b = -0.858, c = -0.136, d = 0.692, 
v = 2.25, and fj, — 3.08. The fitting form is valid for 0.3 < 
m s < 15 keV. This fitting function is shown relative to 
the numerical results in Fig. 01 as well as previous results 
by Ref. 0] • Note that all of the features of the numerical 
results are not obtained in the fit due to the nonthermal 
character of the sterile neutrino distribution, particularly 
for m s = 1.7 keV where peak production occurs near the 
quark-hadron transition. 

The result presented here for the relative sterile neu- 
trino transfer function is similar, yet significantly differ- 
ent from previous work |43j, |4J] , with the difference at- 
tributed the use here of the non-thermal sterile neutrino 
momentum distribution due to the physics described in 
! fTTI The results derived here differ in cosmological pa- 
rameter dependence of T 8 {k) from 2% to 18% and in 
the rapidity of the cutoff fi at 45% relative to that in 
Refs. [43,[44(- Using the transfer function derived here 
and small scale clustering data sets including the inferred 
matter power spectrum from the high-resolution Lyman- 
a forest from Viel et al. Ref. [43 found lower lim- 
its on the mass of the sterile neutrino dark matter at 
1.7 keV (95%CL) from the CMB, the SDSS 3D P g (k) 
of galaxies [63| plus SDSS Lyman-a forest and a 

lower limit of 3.0 keV (95%CL) if the inferred matter 
power spectrum from the high-resolution Lyman-a for- 
est of Ref. jlj is used, which however has significant 
systematic uncertainties. 
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FIG. 4: Shown here are the relative sterile neutrino transfer 
function T s (k) to CDM for the same large scale amplitude 
of perturbations, for the cases of m s = 0.5, 1.7, 8.2 keV with 
increasing wave number k suppression scale, respectively. The 
solid (black) lines are from the full numerical calculation, the 
dashed (red) lines are fitting form in Eq. [11 1 1 , and the dotted 
(blue) lines are the results of Ref. |44|. 



of finite temperature effective mass of active neutrinos 
due to the presence of thermal leptons, and heating of 
the coupled species due to the disappearance of degrees 
of freedom in the plasma. The resulting sterile neutrinos 
have a non-trivial momentum distribution that is grossly 
nonthermal. 

Using the resulting energy distributions of sterile neu- 
trinos, I have calculated the evolution of sterile neutrino 
density perturbations in the early universe through the 
linear regime and provide a fitting function form for the 
transfer function describing the suppression of small scale 
fluctuations for this warm dark matter candidate. The 
results presented here differ significantly from previous 
work due to the inclusion of relevant physical effects. 

The results of the linear perturbation evolution pre- 
sented here are necessary as initial conditions for address- 
ing the questions of structure formation and galaxy for- 
mation in the case of a sterile neutrino dark matter can- 
didate with a non-trivial velocity distribution. Following 
structure formation in this case into the nonlinear regime 
will allow a resolution of the question of whether observa- 
tions of galactic structure at small scales are indicating a 
method of inferring the properties and — ultimately — the 
identity of the dark matter. 



IV. CONCLUSIONS 

Potential problems in galaxy and small-scale structure 
formation indicate the possibility of a small-scale veloc- 
ity damping of perturbations of the type in warm dark 
matter such as sterile neutrino dark matter. Here, I have 
presented the calculation of the production and linear 
perturbation evolution in the early universe of sterile neu- 
trino dark matter. Included were the essential effects of 
the change in Friedmann dynamics dictated by the quark- 
hadron transition and particle population, the alteration 
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